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(^ ■ Abstract 

O 

CN ■ In this paper we present preliminary results from cosmological simulations of modified gravity 

in the dark matter sector. Our results show improvements over standard cold dark matter cosmology. 
The abundance of low-mass haloes in the modified gravity model fit observations better than the 
conventional theory, while the differences of the modified density fluctuation power spectrum differs 
from the standard, ACDM power spectrum are small enough to make these two models observationally 
indistinguishable at large scales. 
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•^z 1. Introduction 

Over the past 30 years cosmology became a successful and mature empirical science. The 
growing accuracy of astronomical observations gave the possibility to test theoretical models. Past 
decade has established the "standard cosmological model" dubbed also as the ACDM cosmology, or 
j->^ . the "concordance model". In this picture the evolution of the expanding universe, dominated by non- 

^D \ relativistic (cold) dark mater particles is governed by Einstein's gravity with a positive cosmological 

t~^ ■ constant. Although the standard model is very successful in explaining many observations it still has 

CD \ some difficulties. The physical nature of the dark matter (DM) sector is still unknown. Moreover, the 

mass distribution in the real universe (in particular, the emptiness of voids) differs significantly from 
the predictions of the ACDM model. Modifications of the standard gravity theory have been proposed 
to deal with this impasse. In this paper we show preliminary results of our work based on one of such 
^ ' modifications. 



2. Long-range scalar interactions in CDM 

The conventional cosmological paradigm is the ACDM model. The model consists of a set of 
assumptions and 'best fit' parameters chosen to describe a certain set of astronomical observations. 
Evolution of the Universe as whole, when we think of scales of order of a few Gpco is governed by 
Einstein's general relativity theory. This evolution obeys Friedman equations. In ACDM the Universe 
is undergoing a phase of accelerated expansion due to large content of Dark Energy, a mysterious 
phenomenon of unknown physical nature that drive global repulsion. On smaller scales (of order 
of hundredths Mpc) we can describe the Universe using Newtonian mechanics. This model is very 
successful on large scales. It agrees well with observations of Cosmic Microwave Background Ra- 
diation (CMBR), dynamics of clusters of galaxies, and provides a natural explanation for the origin 
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and the observed abundances of light elements. Although the power spectrum of the CMBR angu- 
lar fluctuations, predicted by the ACDM model agrees with observations (cf. Wilkinson Microwave 
Anisotropy Probe WMAP results 1 1 1), while the power spectrum of density fluctuations agrees with 
the large-scale clustering seen in galaxy surveys, the ACDM N-body simulations are in contradiction 
with observations on small scales. 

2.1. Shortcomings of ACDM cosmology 

A thorough discussion of all of the difficulties faced by the ACDM is beyond the scope of this 
paper. Instead, we shall briefly point out the most important issues, and in the next section we propose 
a possible solution for them. The problems we address are: 

1. High overabundance of satellites and small galactic haloes in N-body simulations, which 
disagree with Local Group observations, 

2. ACDM simulations suggest that intergalactic and inter-cluster voids are not so quite empty. 
There are still lots of 'debris' - unaccreted matter left in simulated voids. Observations favor an 
opposite picture, where voids are really empty of bright galaxies and matter, 

3. The problem of steep central density distributions of galaxy haloes. Simulations favor a 
cusp with a logarithmic slope for density distribution, whereas observations suggest haloes with a core 
of constant density (see e.g. f2\, f3l), 

4. N-body simulations present a picture of high accretion activity in the recent past. This poses 
a real problem to galaxy formation and evolution theories. Large spiral galaxies would be "dimmed" 
and their spiral structure probably destroyed by late time merging and accretion activity, contrary to 
observations. 

2.2. Scalar „enhancecl" gravity 

Long-range scalar interactions in the DM sector was proposed by P. J. E. Peebles & S. Gubser 
as a possible solution to pure ACDM cosmology problems (cf. ID, Q). Physical motivation for 
this model comes from the string theory. We consider DM particles as strings. Such objects will 
interacts by gravitation and additional long-range force comes from exchange of a massless scalar. A 
screening length for scalar interaction is dynamically generated by the presence of light particles with 
a Yukawa-like coupling to the scalar field. In this paper we consider one species of strings as DM 
particles with a force law rising from the potential of a form: 

^{r) = -G j <fr'^^^{l + Pe^^\ , (1) 

where /? is a "scalar to gravity" ratio parameter and r^ is the screening length. /? is a measure of the 
power of scalar interaction compared to gravity interaction between two DM particles. The screening 
length parameter due to mechanism of dynamical screening is constant in co-moving coordinates. In 
Fourier space this eq. transforms to 
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1 + (fc • rs) 

where (j){k) is the Fourier transformed potential, Hq is present Hubble constant value, Oq is dimension- 
less parameter of present matter density contribution and 5{k) is Fourier transformed density contrast. 
As we can see mechanism of screening leads to screening of scalar force on a particular Fourier modes 
during time evolution. This gives us limiting behavior on large and small scales: 

(j) ^ {I + (3)(t)Newton for r<rs, (3) 

(t>^ <t> Newton for r > r^, (4) 

Considerations on the grounds of the string theory suggest that sensible values of two free parameters 
of this model are /3 ~ 1 and rs ^ IMpc. Earlier numerical and analytical study of this model (see @) 
suggest that presence of the scalar interaction will result in earlier formation of structure on small- 
scales and more rapid pace of formation of structures. Other expected effects are the absence of low 



mass halos in voids and a lower abundance of sub-galactic halos. In the light of presented difficulties 
of the ACDM scenario such properties would be most welcome. 

3. N-body simulations - evolving universe in your computer 

One of the tools of choice for modern theoretical cosmology is a collisionless N-body simulation. 
Computer N-body codes provide a way to investigate non-linear gravitational evolution of complex 
particle systems. There are many variety of N-body codes, algorithms and techniques, but all of them 
essentially model evolution by following the trajectories of particles under their mutual gravity. 

In cosmology one is using N-body codes to model the growth of structures in time from initial, 
small-amplitude inhomogeneities into the highly clustered structure, observed today. The quality of an 
N-body code rests on (1) generating accurate initial conditions (IC) i.e. initial positions and velocities of 
particles, and (2) as precise, as possible solution of the equations of motion to find particle trajectories. 

Using measured power spectrum of CMBR photons we can calculate IC for an N-body experi- 
ment at sufficiently early time of cosmic evolution. Commonly used scheme of generating IC applies 
the Zeldovich approximation to move particles from a Lagrangian point g to a Eulerian point a;(e.g. 
Ha, |18||). 

Having set up the IC we can start to model dynamics of a collisionless system representing 
the Dark Matter particles in the Universe. To do that a N-body code is solving simultaneously the 
collisionless Boltzmann equation (CBE): 

^ + y(y^K.^^ = o] (5) 

dt ^ \ * dxi dxi dvi J 

and the Poisson equation: 

V2$(f) = AirGpir) (6) 

Eq.© is solved using method of characteristics (see |9| for details). The trajectories we get by 
time integrating TV points {fl{t) , v{t)} sampled form the distribution function / at time ij„j form a 
representative sample of / at any time t since / is constant along any given trajectory. The main job 
to be done in N-body code is then to solve Eq. © for a given particle system and move the particles 
to the next time frame according to the equations of motions. 

In this study we are using a Particle Mesh Code (PM) (see IH, 191) that allows one to use a 
N-body simulation code to follow evolution of a system containing millions of particles. 

3.1. The simulation code - AMIGA 

We use a modified AMIGA code (Adaptive Mesh Investigations of Galaxy Assembly), a c code 
for cosmological simulations. AMIGA is a successor of the MLAPM code (Multi-Level Adaptive 
Particle Mesh) (Knebe, Green and Binney (W\) and it is purely grid-based N-body code. Due to large 
number of simulation runs we had to perform, we used a pure PM part of AMIGA code to greatly 
shorten the time of integration at the expense of force resolution. 

To enable AMIGA to do modified gravity simulations, we changed the Green's function form 
to account for investigated effects. The 3-D discrete Poisson equation for standard gravity in Fourier 
space reads 

m = l^^^Hk) X Gfc, (7) 
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where G^, is Green's function: 

G{kx,ky, kz) ^ -tt/ InG"^ sm^ {-Khj^/NG) + sin^ {-nky / N G) + sii?{-Kkz/NG) | . (8) 

Using Eq. Q we modify Green function to get proper potential in scalar-interaction model: 

/^scalar /^Newton , , ( i i P \ /a\ 

We use above modification together with the PM part of the AMIGA code to evolve particle 
system from early epoch to present time of cosmic evolution. Simulations results allow us to study 
the impact of scalar interactions on formation of large-scale structures in the Universe. 



Table 1. Cosmological parameters used in simulation runs 
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4. Simulations results 



To study the impact made by long-range scalar interaction on process of structure formation 
we have performed a series of numerical cosmological simulations. IC for our simulation runs were 
generated using Klypin & Holtzman PMcode |11|. We use adiabatic IC with truncated CDM power 
spectrum generated with a help of cmbf ast code for line-of-sight integration of CMBR (see [.12J ). 
In Table [l] we present values of cosmological parameters used in our simulations. We have tested 
parameter space for [3 = -0.5,0.0,0.2,0.5,0.7,0.8 and rg = 0.5,1.0,2.0,5.0 Mpc for simulations of 
123'^ particles over 256'^ mesh in a box size 200/i^^Mpc. 
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Fig. 1. Power spectrum of DM density perturbations. Top picture shows Pk(k) for a different values 
of [3 parameter and for j^ = l.O/i ^Mpc fixed. Bottom picture illustrates power spectrum for fixed 
value of /3 — 0.5 and for a different values of r^. Lower panels of each picture shows relative ratio 
of measured power spectra to standard Newtonian gravity case. Black dotted line shows the Nyquist 
frequency. Data from simulation of 128'^ DM particles in a box width of 200/i^ Mpc. 



First we analyze the effect of scalar interaction on final power spectrum of DM density pertur- 
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Fig. 2. CMF for redshifts z=0.0 (top figure) & z=1.0 (bottom figure) {(3 parameter is 1.0). Diagrams 
illustrates results from PM simulation of 123'^ DM particles over 16/i ^ Mpc cube length 

bations. FigH] shows the calculated power spectra at redshifti z — from simulations outputs. Each 
curve was obtained by averaging ten simulation runs for a given pair of parameters. For clarity we 
show error bars only for the Newtonian case. On the top picture we keep scalar power parameter /? 
fixed at 0.5 value. This allows us to analyze impact of changing value of screening length (r^). It is 
clear to us that additional scalar interaction in DM substantial changes P{k) only on intermediate range 
of scales related to rg (k ^ l/rs). Increasing the rg parameter increase the amplitude of deviation from 
Newtonian case as illustrated by the lower panel of the top figure. For rg = 5/i~^Mpc case nearly all 
scales in power spectrum are affected. Bottom picture shows P{k) when we keep rg = 1/i^^Mpc and 
change the /3 parameter. We see that changing value of f3 from -0.5 to 1.0 gently increase amplitude 
of fluctuations on scales corresponding to fc ~ 1/iMpc^^. 

Our PM simulations in the box of side length 200/i^^Mpc are limited by force resolution of 
order ^ IMpc and mass resolution ~ 10^"^ Mq. To study the impact of scalar interaction of formation 
of structure on small-scale we have performed a series of simulation of 128^ particles in the box of 
side length 16/i~^Mpc, this corresponds to the force resolution of order ~ SOOkpc and mass resolution 
~ 10^ Mq. In our analysis we included only haloes with more then 20 particles which corresponds 



^ redshift - describes a shift of a light wave towards longer wave tail caused by global expansion of 
the space. It is related to the evolution age of the Universe and the scale factor by a = j^ 



to minimal halo mass M ^ 3.3 x IO'^H^'^Mq. In Figl2] we show calculated Cumulative halo Mass 
Function (CMF) for standard and modified gravity. CMF reflects the abundance of halos and it is 
normalized to number of halos of given mass iV(> M) per cubic megaparsec /i^^Mpc'^. We keep j3 
fixed to 1.0 value and test the Vg value. Top picture is calculated for current time z — the bottom 
picture is calculated at z = 1. The lower panels presents relative ratio of CMF to the Newtonian 
case. These figures show that scalar interaction shifts the merging activity to the higher redshifts and 
towards higher masses. It is also clear that in the presence of scalar interaction with r^ > 1/i^^Mpc 
the abundance of low-mass haloes is substantially lower than in the Newtonian case. We argue that this 
also indirectly implies that inter-galactic voids will be more empty in scalar model due to enhanced 
formation of small-scale structures. To emphasize this, FigU presents redshift evolution of integral 
abundance of haloes with mass M < lO^^h^^MQ. Clearly scalar interaction models show much less 
evolution for 2; > 2 than Newtonian model in this mass range. This indicates that small-scale structure 
growths earlier with scalar-interaction. Undoubtedly in our modified gravity models low-mass haloes 
are absorbed by high-mass cluster and galactic halos more efficiently than in Newtonian case. The fact 
that absorption of low-mass haloes starts on higher redshifts in scalar regime has obvious influence 
on finial lower abundance of this haloes in scalar models. 
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Fig. 3. Redshift evolution of the abundance of haloes with mass M < lO^'^h ^M© 



We turn now to properties of the most massive halo (MMH) in each simulated model. Table [2] 
show general properties of the MMH for tested parameters. We show mass in virialized region, three- 
dimensional dispersion of velocity, velocity at virial radius, the virial radius, triaxiality parameter, 
ellipticities and the spin parameter. We calculate mass, radius and velocity at halo border where 
overdensity reaches 335 value. We define triaxiality parameter and ellipticies as: 



ei 



c 

1--, 
a 



T = 



62 



a 



(10) 
(11) 



where a > b > c are eigenvalues of inertia tensor. For the spin parameter A we use the definition of 
Bullock et al.1131: 

. (12) 
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The most dramatic effect of scalar forces that emerges from parameters in Table [2] is substantial excess 
of 3-dimensional velocity dispersion. The velocity dispersion in scalar models rise more rapidly than 
one could expect from increasing of the halo mass. We account this phenomenon as a purely modified 
gravitation effect. DM particles fall into the steeper potential wells than in Newtonian case and are 
gathered from larger co-moving distances during halo relaxation, thus they obtain more kinetic energy. 



Table 2. General parameters of the most massive gravitationally bound object. 
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As our last but not least result in this paper we want to present calculated MMH's density 
profiles (Figlll) and velocity curves (Fig|5l). We want to outline the striking feature seen in density 
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Fig, 4, Density profile for the most massive halo in scalar and Newtonian simulations. Vertical line 
indicate distance related to the Nyquist frequency of our simulations. 
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Fig. 5. Rotation curves of the MMH. 



profiles, where a trace of solid density core can be found in halos simulated with the scalar interaction. 
Although this tantalizing feature would be most welcomed in context of astronomical observation we 
have to approach it with a caution. In our PM runs we do not resolve reliably scales smaller than 
« 20% of virial radius of the haloes. We have marked distance corresponding to the Nyquist frequency 
and, as can be seen, solid core feature appear below trusted distance threshold. More detailed and 
intense study of this effect is required to decide whether this property is real or just a ghost from 
numerics. 



5. Conclusions 

Our results show that the modification of standard gravity we consider can close the gap between 
some properties of standard ACDM cosmology simulations and today's high resolution astronomical 
observations. Deviations from the standard power spectrum are negligible, and they agree with large 
scale observations. On the other hand, our simulations of scalar interaction show desired clustering 
properties on small scales. While the court on the physical nature of dark matter particles is still out 
it is scientifically justified to probe various theoretical possibilities. In this paper we have shown that 
the long-range scalar DM interaction has the potential to improve the standard cosmological model. 
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